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Abstract 

We address the process of generation of the photon-number entangled states of 
hght in the stimulated nonlinear parametric down conversion process and build the 
simple model describing the generation, not involving the traditional parametric 
approximation. The motion equations for the system of pumping and two-mode 
outgoing field are solved for the case of the strong correlation between two modes, 
the evolution of the state parameters of the generated modes is obtained. The 
solution is briefly analyzed for particular types of photon-number entangled states. 

1 Introduction 

Nonclassical features of the states of light are of utmost interest within the past decades 
as they became the basis for the quantum information theory and it's appHcations, in 
particular quantum computing, teleportation and cryptography All of these are 
based on the fundamental principles of quantum mechanics and mainly involve the es- 
sential resource of quantum entanglement. The recent developments in the field are 
based on the continuous variables information coding |2) into the parameters of multi- 
level quantum systems as distinct from the discrete variables coding with a two-level 
system states. In the optical implementations it mainly corresponds to the consider- 
ation of multiphoton quantum states as the information carriers. At that the special 
attention is given to the entangled multiphoton states, in particular the photon-number 
entangled states of light (PNES), which provide with the strong correlation between 
the intensity measurements of the two spatial modes. Several communication protocols 
were proposed and partially implemented on the basis of the photon-number entangled 
states with the use of intensity difference OJIH or intensity fluctuations coding Q. 
Such states can be generally generated in the nonlinear process of parametric down 
conversion (PDC) l6] E |9l being the well known typical example of the quantum- 
mechanical event when the pumping photons excite the atoms of the nonlinear media, 
which emit photon pairs upon relaxation. While in the first experiments such pairs 
were observed independently, with the increase of process gain in the optical paramet- 
ric amplifiers the generated multiphoton beams started to manifest the nonclassical fea- 
tures, in particular the strong correlation between modes (showing in the sub-shotnoise 
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signal-idler intensity difference fluctuations) and became referred to as the twin beams 
ifTOl [TTl [T2l . At that the modes marginal statistical distributions were thermal. The 
recent advances in the experimental technics (in particular, the use of the powerful 
short laser pulses for pumping and effective multiport photon-counters and CCD cam- 
eras) and usage of cavity-based parametric oscillators enabled several groups to obtain 
strongly correlated states of light which possess Poisson statistics lfT4l[T3l . Meanwhile 
the traditional approach to description of the PDC process is based on the parametric 
approximation, when the pumping mode is considered as classical |7 |. In the present 
paper we build a simple model of the PNES generation without the parametric approx- 
imation, thus keeping nonlinearity in the operator equations. We make the statistical 
assumptions relevant to the considered states and obtain the solutions describing evo- 
lution of the state parameters. At that we assume that pumping is so intense that the 
atoms of the nonlinear media are excited coherently and the PDC process is stimulated. 
This work is the step towards the exact model describing the PDC-based generation of 
the PNES states. 

2 The photon-number entangled states 

The PNES may be written in the Fock number basis as 



where |n, n)) = \n)^^ |n)2, indexes 1, 2 standing for the respective modes and TV is a 
normalization factor. 

Their features are evident thereby: the amplitude mean values for each of the modes 
are zero, while the amplitudes are strongly correlated, i.e. in terms of the modes quan- 
tum operators (ai) ~ (8,2) =0 and (aia2) = A 7^ 0. These statistical properties of the 
states are used within the construction of the PNES generation model. 



Let us describe the possible generation of PNES states by means of a PDC process. 
As it was already mentioned, we assume that the pumping pulses are intense enough 
to coherently excite the whole electron subsystem of the nonlinear media and so the 
pumping energy directly turns to the energy of a generated two-mode electromagnetic 
field excitation. The system can be described by the Hamiltonian: 

H = hiuia'^ai + fujj20'2^'2. + fu-^oO-QCiQ + ihx {ala^aQ — 0-1(120^) (2) 

which includes the interaction term, but doesn't include quantum operators stand- 
ing for the electron subsystem excitation and relaxation and covers only the energy 
swap between the pumping and emitted fields. We assume that the excitation and re- 
laxation takes place immediately, i.e. the model is adiabatic. 

The quantum operators for each of the two modes of the outgoing field satisfy the 
Heisenberg equations: 
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similarly for the hermitian conjugate. The equations lead to ones for the operators 
"1,2 = 2^i,2> which correspond to the photon numbers in each of the two generated 
modes: 

i/l-j-ni = ih^h2 = ihxa2aiaQ + ihxai 0,20,0 (4) 
At that the difference between photon numbers in the two modes vanishes: 

^ {ni - ri2) = (5) 

which means that 

(fii - 112) = (6) 
if initially both of the modes aren't excited. 

Introducing operators A = 0102, A+ = o^Og and iV = n2 + ni we obtain the 
following system of equations: 

= -iujA + X (iV + 1) ao (7) 
^i+=ia;i++x(iV + l)a+ (8) 

j^N = 2x (ia+ + i+ao) (9) 

We suppose that pumping is coherent and the generated state is PNES, i.e. 

(A) = ((^-l A \^)) = A [t); (ao) = {a \a^\ a) ^ a (t); 

(a+) = i+ 1*)) = A (t); (a+) = {a\a+\a) = a (t); 
and that (n'^ = N{X{t)) = N (t). 

Since we are describing the generated state as photon-number entangled, but not 
as the pair of coherent states, we may assume it's statistical independence from the 
pumping state but not the independence between the two modes. The overall system 
state can be presented as \a) (g) j*)), while |^E')) 7^ |q;i) (g) \a2) which is the manifesta- 
tion of modes' entanglement, i.e. (aia2aQ ) = (aia2) (b-q), ^Aa^ ^ = X{t)a (t) and 

Nao^ = N {t)a {t). Thus the system of equations for operators turns to the same for 
the state parameters: 

j^X{t) = -iLuX{t) + x{N{t) + l)a{t) (10) 
p^{t)=iujX{t)+x{Nit) + l)a{t) (11) 

j^N{t)=2x{X{t)a{t) + X{t)a{t)) (12) 

Separating the slow ampUtudes: A {t) = e~'"*A (t); a (t) = e^*'^*a (t) and sup- 
posing that a is real, thus A is real, we obtain the system 
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The solutions for the motion equations ( I13I14I ) has the form 



A{t) = sinhr(i)coshr(i); N = 2sinhV(i), 



(15) 



where 



r{t) = X 



I 



— OO 
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assuming that initial conditions are i^— C)o:T^O,iV^O,A^O. The solu- 
tions describe the connection between the time profiles of the pumping and generated 
modes. The effect of the pumping profile is specified by the integral characteristic ( fT6b . 
in particular for the case of pulse pumping a (<) = {t<0:0;i>r:0;0<<<r: a}: 
T = {t < : 0;t > T : xaT; < t < T : xat}. The result (O is in well concordance 
with the existing description of the PDC process Q. 

4 Particular cases 

In order to establish a link between the model and the physical reality we briefly an- 
alyze it for the two relevant particular types of PNES states, namely the twin beam 
(TWB) and tmo-mode coherently-correlated (TMC) states. 

The twin beam states which are also referred to as squeezed vacuum states can be 
presented in the form ([T]) as 



where generally x G C and < < 1, but without loss of generality we assume x 
as real . TWB states are Gaussian and are widely used for CV teleportation and dense 
coding. The marginal statistical distributions of the TWB modes are equal to thermal 
states. 

The TMC states being the degenerate pair-coherent states ifTSl can be given as 



where Ia{x) is the modified Bessel function. These non-Gaussian states are eigen 
for the product of the modes' annihilation operators 0102 |A)) — A|A)) and possess the 
sub-Poisson photon-number statistics for each of the modes. 

In order to estimate which of the states is closely to be described by the model, 
while the mode amplitudes have zero mean values we introduce the observables which 
we call pair-quadratures similarly to the mode quadratures and construct them of op- 
erators A ~ aia2, = (1^0,2 as (7+ = A + A+ and C- = We estimate 



the evolution of a pair-quadrature dispersion Dc — (&) ~ (C) , where C = C± 
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within the model directly diffrentiating the dispersion obatined with averaging by the 
time-dependent solutions: 



dPc 
dt 



odel 



-(((vI/(t)|C7^|M^(t)))-((vI/(t)|C|*(t)))- 
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assuming that the state evolution is decribed by the time dependence of state pa- 
rameters ( fTSl l. We compare this to the exact evolution of the pair-quadrature dispersion 
obtained directly by averaging the pair-quadrature derivatives by the particular states 
either TWB ^ or TMC dH: 



dPc 
dt 
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where ^ is obtained from the Heisenderg equation ^ 



[C, H] for the system 
Hamiltonian 

For the C+ pair-quadrature the resulting derivative within the "exact" approach is 
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where Q = oo + oo is the pumping quadrature. Calculating this derivative for the 
particular type of states (TWB or TMC) and comparing it to the expression obtained 
by calculation through averaging by time-dependent solutions within the model, we get 
that for TWB 



dD, 



dt 



TWB 



dDc+ 
dt 



TWB 



odel 



8xa 



x(x^ + 1) 



(l-a;2)2 
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while for TMC 



dD, 



dt 



TMC 
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dt 
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4xAa, 
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thus the built model is more likely to describe the TWB states generation. 



5 Conclusions 

In this work we built the simple model describing the photon-number entangled states 
generation in the stimulated parametric down conversion process. At that we didn't 
make the traditional parametric approximation and considered pumping mode as quan- 
tum. The special statistical properties of the states enabled obtaining the solution for 
the motion equations describing the process of generation in terms of the evolution 
of the sate parameters. The solution represents the relation between time-profiles of 
the pumping and generated modes. The model is briefly analyzed with respect to the 
particular types of photon-number entangled states and is more likely to describe the 
generation of a twin-beam state. 
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